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On Light Elliptically Polarised by Pejlexion, especially near the 
Polarising Angle : a Comparison with Theory.* 

By Richard C. Maclaurin, M.A., LL.D., late Fellow of St. John's College, 
Cambridge, Professor of Mathematics, Wellington, New Zealand. 

(Communicated by Professor J. Larmor, Sec. R.S. Received January 14, — Read 

May 23, 1907.) 

Taking the phenomena of light to be due to periodic displacements in a 
rotational ether, we shall investigate the exact character of these displacements 
in the case of reflection from an isotropic transparent substance. Two well- 
known vectors play the principal parts in any such optical discussion — the 
displacement of the medium and its curl ; in the electro-magnetic theory the 
former is proportional to the magnetic force and the latter to the electric 
" displacement/' To give precision to the discussion we shall deal throughout 
with the displacement of the medium, there being no difficulty in dealing in 
a similar fashion with the curl. 

Let the displacement in the incident light be of unit amplitude and 
polarised in a plane making an angle 6 with the plane of incidence. It can 
be resolved into two components £ = cos 6 cospt in the plane of incidence, 
and 7] = sin 6 cospt perpendicular to this plane. The corresponding com- 
ponents in the reflected beam are then given by the formulae 

£= — -R' cos#eos(^— //) and i) = — R sin 6 cos (pt — p). 

* [With reference to an oversight in his paper " On Metallic Reflexion," ' Roy. Soc. 
Proc.,' vol. 77, Jan., 1906, which has recently been pointed out privately by Mr. C. E. 
Hasse, Prof. Maclaurin writes that he had himself observed it immediately on receipt of 
the printed paper, but that fortunately the blemish does not invalidate any of the results, 
beyond requiring a slight modification of the argument in one place. 

At the foot of p. 218, in quoting the formula E = p~ 2 \ \xdx from the previous paper, 
the first factor /x~ 2 was left out. When it is inserted, the modulus of E is no longer 
large, and the following modification becomes necessary. We have 

A,4-a, A, - 2cd, cos <b f , n N m ^ ^~ 

: a* T / = T 1 ' H"5 7a7—TZ' {(E- 1W 2 + (F-E) sin 2 c&} ; 

A. 1 + a l Ai L f* 2 cos 2 9 — cos 2 9 lv ; r v ' J 



A, 



r-2 



-O-O 



■ 1+ _2^cosj$_ (E _ 1} ,-|. 

IX 1 cos 2 <p — cos 2 q> J 



/A 2 cos 2 <p f — COS 2 (j> 

If both the terms (E — 1) ix 2 and (F — E) sin 2 cj> are to be retained, we would require two 
complex constants to specify the layer of transition, which the experimental knowledge 
would hardly be adequate to determine. When, however, (f> is small, the term (F — E) sin 2 <f> 
Avill be negligible in comparison with (E — 1) /a 2 ; and even if <jf> is not small, the first term 
will not usually be more than one-tenth of the second, the modulus of the factor fi 2 in the 
latter being large. The neglect of the first term may therefore be expected to produce an 
error of about one-tenth in the estimated correction due to the layer of transition for large 
incidences, and practically no error at all for small incidences. This degree of approxima- 
tion fortunately remains adequate, in relation to the exactness of the data utilised in the 
argument.] 

VOL. LXXIX. — A. 2 L 



482 Prof. R. C. Maclaurin. [Jan. 14, 

Here R' and E are given by the sine and tangent formulae of Fresnel, viz. : — 

R , _ sin (eft -eft') . E _ + tan(j>-^ / ) 
sin (0 + <//) ' ~ tan (<j> + <j>') ' 

where <j> and <£' are the angles of incidence and refraction, and the upper or 
lower sign is chosen so as to make R positive throughout, p and p are the 
changes of phase produced by reflection. In the case of an abrupt transition 
from one medium to the other, p f would be zero throughout, while p would be 
zero or ir according as <£ were less or greater than Brewster's angle tan -1 /*. 
In the more usual circumstances of a gradual transition p and p are given 
very approximately by the formulae : — * 

tan p- ch x V C ° S * [(F ~ E) Sin2 *,- (1 ~ E) l3 , 

^ (/jb 2 — 1) [sin 2 ^>—/>6 2 cos 2 <5f)] 



£6 J — 1 



and the difference of phase A = p — p by the formula 

tanA = rflX 2/i'{(F-l)-^(l-E)} sinjjan^ 

yL6^ — 1 tan J <£> — yLt J 

_ a p? sin ^> tan <£ 

~ ~ ° (l + /x 2 )* ' lan 2 ^^ 2 • 

Here d\ = 2'ird/\ i where X is the wave-length in air and d is the thickness 
of the transition layer; E and F are numbers depending on the law of 
distribution of the refractive index within the layer, and e is the coefficient 
of ellipticity at Brewster's angle. 

The formula for tan p f shows that it is proportional to cos <£, so that it will 
diminish to zero as the incidence increases. In nearly all cases that have 
been examined E is less than unity and d\ is a very small quantity. Thus p 
will be small and positive throughout, approaching the limit zero as the 
incidence becomes grazing. At direct incidence, where there is no distinc- 
tion between the plane of incidence and that at right angles to it, p is 
equal to p . 

With the substances that are dealt with specially in this paper, and in 
nearly all cases where the experimental determinations are well assured, p 
increases with the incidence, slowly everywhere except near Brewster's angle, 
in the neighbourhood of which its increase is very rapid. In such circum- 
stances p is greater than p and the crests of the waves representing £ are 
reached sooner than the crests of the waves representing 77. Thus tj, the 

* See ' Koy. Soc. Proc., 5 A, vol. 76, 1905, pp. 55 and 57. 
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displacement perpendicular to the plane of incidence, lags behind £, the 
displacement parallel to that plane, the difference of phase being A.* 

MS/e have seen that the components of the displacement in the reflected 
beam are given by 

£ = — K' cos 8 cos (pt ~p) and tj = — E sin 8 cos (pt—p). 

If = or -|7r, then r\ = or f = 0, so that the reflected light is plane 
polarised in or perpendicular to the plane of incidence. In other cases, 
in general, it is elliptically polarised, the elliptic orbit being obtained by 
eliminating t from the above equations for f and tj. In this way we obtain 
as the equation of the ellipse 

& t] 2 2£t; cos A 



+ zzTr-r-z-z — - ^^, . — 7i ^ = Sir A. 



E' 2 cos 2 8 E 2 sin 2 8 EE' sin 8 cos 

In the special cases of direct or grazing incidence the ellipse degenerates 
into a straight line. If A = the reflected light is plane polarised in an 

E 

azimuth 8' given by tan 6 f = -? = ^, tan # = e tan #, while if A = ir the 

light is plane polarised in an azimuth 6', where tan 9' = — etan#. These 
two cases would represent the state of affairs below and above Brewster's 
angle respectively, if the transition from one medium to the other were 
abrupt, so that in such circumstances the reflected light would be plane 
polarised throughout. 

By considering the signs of £, tj, d%/dt and drj/dt in the formulae above, 
it appears that the motion in the elliptic orbit takes place in a counter 
clockwise sense. To give precision to this statement we may express the 
matter thus : — Look along the ray (whether incident or reflected) in the 
direction of the propagation of the light. Turn counter clockwise from 
the plane of incidence through an acute angle 8 until the direction of the 
incident vibrations is reached. The motion in the orbit is then clockwise.f 

"We proceed to consider the elements of the ellipse referred to above. 
Putting tan a = e tan 8, the equation of the ellipse becomes 

£ 2 tan 2 u + t] 2 — 2f?7 tan a . cos A = E 2 sin 2 A sin 2 8. 

If ty be the azimuth of the major axis, we have 

, . 2hf 2 tan « cos A , n A 

tan 2-vir = — — -. = _- — - — .- — - = tan 2u . cos A. 

a —b 1 — tan J a 

* Of course, if p' were greater than p, this statement would have to be reversed. 
Moreover, if we were dealing with the curl instead of the displacement, (Rp) would be 
interchanged with (K'/Z), so that the component parallel to the plane of incidence would 
lag behind the other. 

t If we were dealing with the curl, instead of the displacement, the motion would be 
counter clockwise. 

2 L 2 
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The eccentricity (e) is given by 

e*> = 1 4:(a'b'-h' 2 ) 

(2-e 2 ) 2 {a' + Uf 

Putting sin 2^ = sin A sin 2 a, this gives 



1— sin 2 2a sin 2 A. 



e = \Acos_2t) 
cos 7 



while for the semi-axes (<x and &) we have 

E sin cos r/ , 

a = ; '- and 



sin« 



a 



= tan 7. 



We shall develop various laws regarding the variation of the elliptic 
elements from these simple formulae, and by way of numerical illustration 
shall apply them to the case of reflexion from diamond and realgar. These 
substances are selected, because, for them, we have the careful experiments 
of Jamin to compare with theory, and so to assure us of the firm foundations 
on which w r e are building, and also because while their refractive indices 
are almost identical, their coefficients of ellipticity at the Principal Incidence 
differ very considerably, so that their comparison should bring into strong 
relief the influence of this ellipticity. The ellipticity in the case of realgar 
is abnormally high ; it stands near the head of Jamin's list, while diamond 
occupies a place near the middle. In realgar we shall take fi = 2*420 and 
6 = 0*1120, and for diamond jjl — 2*434 and e = 0*0242. This makes 
0o = tan -1 yu, = 67° 33' in the case of realgar, and 67° 40' in the case of 
diamond. The difference of phase (A) and the coefficient of ellipticity 
(e = K/B') are then as follows : — 



i $ = 0. 

i 


10°. 


30°. 


60°. 


65°. 


67° 30'. 


67° 40'. 







0°18' 


3° 0' 


27° 43' 


57° 6' 


89° 18' 


91° 44' 


1 

A (diamond) ... 





0° 4' 


0°39 / 


5° 54' 


17° 18' 


77° 26' 


90° 




1 


-9750 


-7829 


'2022 


'1343 


'1120 


-1125 


e (diamond) ... 


1 


-9750 


7838 


'2052 


'0749 


'0247 


-0242 



■ 


<p = 0. 


68°. 


70°. 


75°. 


80°. 


90°. 







96° 43' 


123° 14' 


155° 56' 


168° 


180° 


A (diamond) ... 





111° 24' 


160° 24' 


174° 6' 


177° 21' 


180° 




1 


'1140 


-1352 


-2370 


'4329 


1 


€ (diamond) ... 


1 


0-0281 


-0699 


-2348 


'4301 


1 
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These results are represented graphically in figs. 1, 2, 3, and 4 below. In 
these figures the crosses ( x ) indicate the results of Jamin's experiments, so 
that the agreement between theory and experiment is readily estimated. 
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Fig. 1. — Difference of phase for realgar. 
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Fig. 2. — Coefficient of ellipticity for realgar. 
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We have seen that the azimuth of the major axis of the ellipse is given by 
the formula tan 2-^ = tan 2a . cos A, where tan a — e tan 6. 

At direct incidence e = 1 and A = 0, so that i|r == a = 6 ; at grazing 
incidence e = 1 and A = ir, so that i/r == — a = — 0. Thus as the incidence 
increases ^ will diminish from 6 to-— 0, and for most azimuths it will pass 
through the value zero at the Principal Incidence where cos A = 0. 
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Fig. 3. — Difference of phase for diamond. 
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Fig. 4. — Coefficient of ellipticity for diamond. 



However, for azimuths greater than XI = 90°— wo, where e = tan o>o, 
this will not be the case. In such circumstances tan yjr will pass through 
infinity at the Principal Incidence, and ^ will change from 90° to — 90°. 
Since tan 2i|r = tan 2a . cos A, it will be seen that if <£ be less than the 
Principal Incidence, tan 2y}r is negative, and yfr is less than 90° ; while if <£ be 
greater than the Principal Incidence, i|r is greater than— 90°. Hence, for 
such azimuths the graph for -^ will, in the neighbourhood of the Principal 
Incidence, have the form represented in figs. 5 and 6 below. 

There must, therefore, be some incidences for which the azimuth yjr is 
a minimum or maximum. To find them we note that ijr is stationary when 
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tan A . ~ = tan6> ~- . ~. The minima will be found in the neigh- 

dcf) tan a . cos 2a, d<j> 

bourhood of the Principal Incidence (<£ ). Putting <j> = <£ 4-^, we have 

tanA = A sin f tan t = ~(l-h^ + ^ 2 +...); 
tan J <£ — /r a; 

therefore tan A ^ = ( / V^ 7 \ . 

rf^ a? (1 4- 26a?) 

Also € = eo(l— 2psx+px 2 ). This makes e a minimum when a; = s, and, 
as the minimum is very near O > s will be very small. 

il=2pe (-s+aO; 

cos 2* = 1 "- tai12 ^ = l-tan 2 ao(l~-2^^+^ 2 ) 2 _ cog 2ao (1 + ^ + q X 2 \ 
1 4- tan 2 a 1 4- tan 2 a (1 — 2psa? +^ 2 J 2 

where 5' is a small quantity. 

Thus the equation for x becomes 



a; (1 4- 26a?) cos 2a (1 -4- s'a; + qx 2 ) (1 — 2psa; +_pa? 3 ) * 

p is positive, so that for real (small) values of x we must have cos2a 
negative. This requires « >45° and tan > eo" 1 or > 90° — coq >X2. 

It is then only for such azimuths that yjr is stationary anywhere. 

From the equation tan a = e tan we see that « increases from 45° to 
90° as increases from O to 90°. Hence in this range cos 2« varies from 
to — 1, and we may put cos 2« = — 2pK 2 , where k 2 is small when 6 is 
near O. The equation for x then becomes 

x 2 — sx = a; 2 (1 — feu) (1 -f s'x 4- £a? 2 ) (1 — 2psa? +p% 2 ). (1) 

If a;, s, and /c be all of the same order of small quantities, then by retaining 
only the lowest terms in (1) we get x 2 —sx = # 2 , so that x = ^s±^/(/c 2 + ls 2 ). 
One root is positive and the other negative, the positive root being the 
larger. Thus the maximum beyond the Principal Incidence is somewhat 
further from that angle than the minimum that occurs before the Principal 
Incidence is reached. 

If s be so small as to be negligible, then s = s' = and (1) becomes 
x 2 = # 2 (1 — bx) on neglecting terms of the fourth order of small quantities. 

This gives a positive root, x = #(1 —■£&#), and a negative root, 
x = -/c(l4^4 1^ & be negative, the positive root is numerically the 
larger ; while if b be positive, the negative root is the larger. It will be 
found that b = (3— yu, 2 )/2yu, so that the sign of b depends on whether /jl 
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is greater or less than ^3. It will be negative for such highly refringent 
substances as diamond and realgar, but more usually positive. 

If s be not* zero, but small compared with x, we get for the positive root of 
(1) above, x? = k 2 (l — bfc) + s/c, so that x = k + ^s — |&/c 2 , and for the negative 
root x 2 — k 2 (l + b/c) — s/c and x = — k-{-\s — j^b/c 2 . 

The following tables give the values of ty for the azimuth 6 = 89° in the 
neighbourhood of the Principal Incidence for diamond and realgar. It will be 
seen that the maxima and minima are much nearer the Principal Incidence 
in the first case than in the second. The figures that follow represent the 
same results graphically, fig. 5 referring to diamond and fig. 6 to realgar ; the 
scale of in the former figure being four times that in the latter : — 

Diamond. 



xb 


65°. 


67° 30'. 


67° 35'. 


67° 39'. 


67° 40'. 


77° 24' 


74° 15' 


72° 38' 


86° 6' 


90° 




-90° 



4>. 

■ 


67° 41'. 


67° 45'. 




•J/ 


: -86° 17' 


-70° 54' • 









67° 50'. 


68°. 


70°. 


-64° 21' 


-69° 21' 


-76° 41' 



Eealgar. 




cj>. 


68°. 


69°. 


70°. 


71°. 


72°. 


73°. 


74°. 


75°. 


* 


-88° 57' 


-87° 2' 


-85° 54' 


-85° 28' 


-85° 24' 


1 
-85° 34' 


-85° 52' 


-86° 9 



From these tables and figures it appears that, in the case of diamond, the 
maximum value of \(r occurs about 10 minutes beyond the Principal Incidence 
and the minimum about 5 minutes below that angle ; while the corresponding 
quantities for realgar are 5 degrees and 4-J degrees respectively. 

Eor azimuths less than O we have seen that there are no such maxima and 
minima, and that i/r decreases throughout from 6 to — 6 as the incidence 
increases, and passes through zero at the Principal Incidence. The following 
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Fig. 5. — Diamond, 
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Fig-. 6. — .Realgar. 



490 



Prof. It. C. Maclaurin. 



[Jan. 14. 



table sets out the values of ^ for various azimuths (6) and incidences (<£) in 
the case of diamond and realgar. The first entry in each case refers to realgar 
and the second to diamond. 



0. 


<p = 10°. 


30°. 


60°. 


65°. 


67° 30'. 


67° 40'. 


68°. 


70°. 


75°. 


80°. 


o 


o / 


o / 


o / 


o / 


o / 


o / 


o / 


o / 


o / 


o / 


K, { 


9 45 

9 45 


7 52 
7 52 


1 49 

2 3 


44 
44 


1 
3 


-02 



-08 
-06 


- 45 

- 40 


- 2 12 

- 2 21 


- 4 16 

- 4 20 


30 | 


29 22 
29 23 


24 19 

24 21 


5 54 

6 45 


2 25 
2 22 


3 
11 


-07 



- 27 

- 19 


- 2 28 

- 2 11 


-78 
- 7 41 


-13 46 
-13 56 


45 j 


44 16 

44 16 


38 4 
38 5 


10 18 

11 30 


4 14 
4 6 


5 
19 


- 12 



- 46 

- 33 


- 4 18 

- 3 46 


-12 19 
-13 9 


-23 6 
-23 15 


60 | 


59 22 
59 22 


53 36 
53 37 


17 38 
19 30 


7 29 
7 5 


8 
32 


- 21 



- 1 23 

- 57 


- 7 36 

- 6 31 


-21 1 
-22 3 


-36 42 
-36 40 I 


n { 


74 38 
74 17 


71 6 

71 8 


36 3 

37 25 


18 
15 2 


21 

1 10 


- 53 



- 3 28 
-23 


-18 27 
-13 32 


-41 10 
-41 13 


-58 27 
-58 6 


85 j 


84 53 
79 45 


83 37 

77 19 


68 21 
67 


64 33 
40 22 


88 35 
3 48 


-86 36 



-78 11 
- 6 45 


-64 41 
-38 14 


-71 
-69 42 


-78 48 
-78 31 


89 j 


88 58 
88 59 


88 44 
88 43 


85 38 
85 9 


85 54 

77 24 


89 53 

74 15 


-89 43 
90 


-88 57 
-69 21 


-85 54 
-76 41 


-86 9 
-85 46 


-87 44 
-87 40 



These results are exhibited in figs. 7 and 8 below, the first dealing with 
realgar and the second with diamond. 




Fig. 7. — Eealgar. 
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Fig. 8. — Diamond. 

We have seen that if the transition from one medium to another were 
abrupt instead of gradual, the reflected light would always be plane polarised. 
The azimuth of the plane of polarisation would be a or — a, according as the 
angle of incidence were less or greater than the Principal Incidence. The 
angle a is given by the equation tan a = € . tan 0, and it must be remembered 
that € is not the same in the neighbourhood of the Principal Incidence in the 
two cases of abrupt and gradual transition. Fig. 9 represents an ellipse whose 




■*! 



Fig. 9. 



centre is O and of which OA is the semi-major axis. OI and- OP are parallel 
and perpendicular respectively to the plane of incidence. EH and FH are 
tangents parallel to OP and OI respectively. OE represents the amplitude of 
the displacement parallel to the plane of incidence, and OF the corresponding 
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amplitude for the perpendicular plane. Thus OE = E / cos 6 and OF = E sin 8, 
tan HOE = OF/OE = etan 9, so that HOE represents the angle a. It will 
be seen from such a figure, or from the equation tan 2yfr = cos A x tan 2a, that 
ty is numerically less or greater than a, according as a is less or greater than 
45°. Let % represent the amount of turning towards the plane of incidence 
in passing from the azimuth of plane polarised light in the case of abrupt 
transition to the axis of the elliptic orbit when the layer of transition is taken 
into account. % will be positive or negative according as a is less or 
greater than 45°, i.e., according as e tan 6 is less or greater than unity. Thus 
X will be negative for all angles of incidence when the azimuth is greater 
than XI. It will also be negative for high and low incidences, when e is 
nearly unity, for azimuths greater than 45° ; in all other cases it will be 
positive. At the Principal Incidence, % is zero for azimuths less than X2, 
and 90° for greater azimuths. Thus the layer of transition not only replaces 
linear by elliptical polarisation, but, when the azimuth exceeds X2, it turns the 
direction of the " major axis " from the plane of incidence through a right 
angle.* It will be seen later that the polarisation is circular at the azimuth 
O, so that there is no discontinuity in passing through that angle. The 
following table sets out the values of % for various azimuths and incidences, 
the first entry referring to realgar and the second to diamond : — 



e. 



= 10°. 



30°. 



60 c 



65°. ! 67° SO 7 . ! 67° 40' 



68°. 



70 c 



75°. 



80 c 



10 



{ 



SO- 



Ol 



45 



60 



75 



85 



89 




















































12 





1 


0' 


45 





2 


1 


12 


o- 


6 


1 


40 





5 





58 





'2 


-1 


45 





5 


-0 


33 


1 ■ 


1 









0*2 





0-2 


o- 


0-7 


o- 


0-8 


o- 


0-11 


o- 


0-25 


o- 


0'19 


o- 


0-23 


o- 


0*35 


o- 


25-17 


-87- 


0-12 


o- 


9*36 


-85- 


0-31 


-58- 



1 
1 

1 

2 

1 

2 

5 
6 

47- 
34 

50 
26 












o- 








0-2 


o- 


o-i 


o-i 


0-3 


o- 





o-i 


0'8 


o- 





o-i 


0-7 


1- 





o-i 


0-14 


o- 


0-2 


0*4 


0-23 


1- 





0-3 


0*23 


o- 


0*11 


0-38 


1-12 


o- 





0-5 


1-5 


o- 


-84-26 


-69-34 


-24-48 


-1- 





0-15 


0-23 


o- 


-79-1 


-51*42 


- 9-24 


-o- 


-90-0 


-39-36 


- 0-42 


-o- 



•12 
•1 

•39 
•2 

•1 

•4 
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•5 
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•1 

•16 

•8 

•22 
•1 
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0-16 

o-i 









•0 


•0 




19 
1 

10 
1 

13 

1 

13 
1 



-0*3 




We shall next consider the variation of the quantity bfa, the ratio of the 
axes, on which the eccentricity of the ellipse depends. We have b/a = tan y, 

* See figs. 21 and 22 below. 
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where sin 27 = sin A sin 2a and tan a = e . tan #. Thus b/a vanishes when 
= or -J-tt, or when A = or tt. These, of course, are the cases of plane 
polarisation, 6 = or |-tt being azimuths in or perpendicular to the plane of 
incidence, and A = or it being direct or grazing incidence at any azimuth. 

For any given azimuth b/a begins at zero for direct incidence, and ends at 
zero for grazing incidence. It must pass through certain maxima or minima 
in this interval. From the equation for 7 we see that b/a is stationary when 
cos A . dA sin 2 a 4 2 sin A cos 2a . da = 0. 

This is satisfied at the Principal Incidence where cos A = 0, and da = 0, 
since a is stationary with e. As we move in either direction from the 
Principal Incidence, sin A diminishes and, for all azimuths less than XI, 
sin 2a increases. For substances such as diamond, for which the variation of A 
is very rapid in the neighbourhood of the Principal Incidence, the variation 
of 7 will be mainly controlled by that of A, and we shall have b/a a maximum 
at the Principal Incidence. With other substances, such as realgar, for which 
the variation of A is less rapid, the increase in sin 2a may more than counter- 
balance the decrease in sin A. In such circumstances b/a will be a minimum 
at the Principal Incidence, and there will be two maxima on each side of 
this and not far from it. Especially with substances of high refractive index, 
the increase of e beyond the Principal Incidence is much more rapid than its 
decrease before, so that the maximum beyond the Principal Incidence will be 
further from it than the one below that incidence. For azimuths greater 
than O, we have 2a > 90°, so that sin 2a diminishes with the increase of $. In 
this case, therefore, b/a diminishes as the incidence increases beyond the 
Principal Incidence. There can thus be no maximum beyond that incidence, 
the only maximum possible being one very near the Principal Incidence and 
below it. 

The following tables set out the values of b/a and of the eccentricity (e) 
for realgar and diamond, and the results for b/a are exhibited graphically in 
figs. 10 and 11 : — 
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Fig. 10. — Eealgar. 
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Fig. 11. — Diamond. 



On p. 493 are considered the variation of b/a for different incidences 
at a given azimuth. Its variation for a given incidence at different azimuths 
is very simply discussed. Under such circumstances e and A are given, so 
that sin 2j varies as sin 2a, and 7 is greatest when a = 45°, in which case 
7 = -|A or -J (tt— A), according as the incidence is below or beyond the 
Principal Incidence. Thus at a given incidence hj a begins at zero with 
6 = 0, and ends with zero at 6 = 90°. It passes through a maximum at the 
azimuth 6 = cot~ 1 e J and its maximum value is tan JA or cot -JA, according 
as the incidence is less or greater than the Principal Incidence. The 
following table gives the position and magnitude of the maxima for various 
incidences : — 
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4>. 


10°. 


30°. 


60°. 


65°. 


67° 30'. 


9 (realgar) ......... 


45° 44' 


51° 56' 


78° 30' 


82° 21' 


83° 37' 


hja (realgar) 


-00262 


-02619 


-24686 


-54409 


-98786 




45° 42' 


51° 55' 


78° 24' 


85° 43' 


88° 35' 


hja (diamond) ... 


0*00058 


-00567 


-05153 


-15213 


-80163 



*. 


67° 40'. 


68°. 


70°. 


75°. 


80°. 


9 (realgar) 


83° 35' 


83° 30' 


82° 18' 


76° 40' 


66° 35' 


hi a (realgar) 


-97020 


-88888 


-54032 


-21316 


-10510 


6 (diamond) 


88° 37' 


88° 30' 


85° 


76° 47' 


66° 43' 


hi a (diamond) ... 


1 


-68215 


-17273 


-05153 


-02328 



These results are represented in figs. 12 and 13 below, the first repre- 
senting the plan (0) and the second the elevation (b/a). 




—~p 10 



Fig. 12. — The continuous curve represents the results for realgar, the dotted one for 

diamond. 



VOL, LXXIX. — A. 
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Fig. 13. — The continuous curve corresponds to realgar, the dotted curve to diamond. 

The axes of the elliptic orbit are given by the equations 

_ E sin cos 7 __ B/ cos cos 7 
sin a cos a 



and 



, , E sin sin 7 E' cos sin 7 
= a tan 7 = : f - = ; * 



sm a 



cos a 



At direct or grazing incidence we have 7=0 and a = 0, so that 
<x = E = E' for all azimuths. We shall consider how a varies with the 
incidence and azimuth. We have 

ah E 2 sin A sin 20 



a 2 = 



tan 7 



2e tan 7 



Thus for a given incidence a 2 = a: . -7— — , where /e is a constant. 



we have 



2a da 



K 



tan 7 



= cot 2 7 [2 tan 7 cos 20 d0-~ sec 2 7 sin 20 $y] 



Hence 



= cot 7 . d0 



cos 20 



cos 2 



a 



cos 27. 



Since tan a = e . tan 0, we have a< and cos 2a > cos 20, 2a and 20 both 
being less than tt. Moreover, 2y is never greater than ^7r, so that cos 27 is 
positive and less than unity. Hence if 2a be less than ^7r, so that cos 2a is 

•u COS 2a ^ ozl oz j COS 2a . J ,, 

positive, we have -~->eos2a > cos 20 : cos 20 — is consequently 

* ' cos 27 cos 27 u J 

negative, and a diminishes as increases. However, this argument is no 

longer good when 2a is greater than 90°. We shall, therefore, proceed to 



1907.] On Light Elliptically Polarised by Reflexion. 499 

investigate the roots of the equation cos 20 = cos 2a/cos 27. Combining 
this with sin 2y = sin A . sin 2a, we get — 

cos 2a 



cos 20 



= cos 27 = +(1 —sin 2 A sin 2 2a)*, 



the positive sign being taken, as cos 27 is positive. For brevity, write 
x = tan 2 a, 5 = sin 2 A, c = e 2 , so that we have 

(1— a?)(c + a?) __ /., 4«e \* 
(l+a?)(c-a?)"" \ (1+^)7 

an equation obviously satisfied by # = 0, corresponding to a = = 0. 
Squaring, and clearing of fractions and of the root x = 0, we get 

x*(l — c— s) + 2scx— <? + £ 2 (l— 5) = 0; 

x(l — c— s) = — sc±(l— c) {c(l— s)}*. 

The process of squaring introduces an irrelevant root. To exclude this 

the sign of . — ~ is to be positive, and as — is necessarily positive 

± ~j- X C—~X 1. -{- X 

this requires to be positive. The value of x given above makes 

<?— x ° 

I—- a? __ ~- /l— s\*__ ~ cos A 

C — X \ C ) € 

and as cos A is positive or negative according as the incidence is less or 
greater than the Principal Incidence, we must take the lower sign below the 
Principal Incidence and the upper sign beyond it. 
Thus below the Principal Incidence we have 

tan 2 6 — ^ tl2 ^ — x — •_ 1 — 6 cos A 
e 2 e 2 e(cosA — e)* 

But 1 — 6 cos A is positive, and cos A — e is also positive, except very near the- 
Principal Incidence. Hence, except in the immediate neighbourhood of the 
Principal Incidence, tan 2 would be negative if a were stationary, so that 
there is no real value of for which a is stationary. Beyond the Principal 
Incidence our formula gives 

tan»fl = 1 / +{CMA . 
e (cos A -f- e) 

The numerator of this fraction is positive and the denominator negative, 
except very near the Principal Incidence. We see then that da is negative 
throughout, except in the immediate vicinity of the Principal Incidence, so 
that a diminishes as the azimuth increases for all incidences outside a small 
region near the Principal Incidence. This region extends on each side of the 
Principal Incidence from cos A = e to cos A = — e. Suppose that 

2 M 2 
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cos A = ±e when <j> = <^>o+^> where <£o is the Principal Incidence and x is 
small. We have 

tan A = -4e„ -^^ pt£±± = _ 2e „ (-g-V (j + 1±^T 

^(l+yct 2 ) tan 2 <£— /z 2 \1+A* 2 / \ /* / 

so that 

cos A = cot A (approx.) = — ~ I ( ^ - j = + e = + eo (1 +px I qx 2 ), 

where, as we have seen, p is very small. 
This equation gives, very approximately 

x= ± 2e ^ ^ 



l + fi 2 /' 

and owing to the smallness of e 2 this is usually very small. For diamond 
we have x = + 2' 55" and for realgar x = + 1° 2' 55". 

At the Principal Incidence a is stationary when x = 1, ix. t when 
cot = eo and # = X2, the azimuth giving circular polarisation. The formulae 
for a and b are 

__ E sin # cos 7 __ B' cos cos 7 ^ 7 __ E sin sin 7 __ B 7 cos sin 7 
sin a cos a sin ot cos a 

where sin 2 7 = sin A sin 2 a. 

At the Principal Incidence sin A = 1, so that 7 = a or \tt— a. If 7 = a 
we have a = B' cos and & = B sin#, and we must interchange a and h if 
we take 7 = J77-— a. As <x is the greater of the two semi-axes, the first or 
second of these values must be assigned to 7, according as B 7 cos is greater 
or less than B sin 0, i.e., according as is less or greater than X2. 

Thus, at the Principal Incidence, if be less than X2 we have a = B' cos 0, 
while if be greater than X2 we have a = B sin 0. Hence, after passing the 
azimuth X2, a increases from Bo (1 + eo 2 )"^ ^° ^o- 

It is, therefore, a minimum at the azimuth XI. It must be remembered, 
however, that in passing through this azim uth, the angle y}r which the major 
axis makes with the plane of incidence suddenly changes from 0° to 90°. 
This means, of course, that the longer axis is perpendicular to the plane of 
incidence, instead of being in that plane. At the Principal Incidence the 
axis in the plane of incidence is Bo 7 cos throughout, while the axis 
perpendicular to that plane is B sin 0. As the azimuth increases, the first 
axis diminishes steadily from Bo 7 to zero, while the second increases steadily 
from zero to B . 

We have still to consider the variations of a with the incidence for a given 
azimuth. Prom the formula for a, viz., a = B sin cos 7/ sin u = B' cos cos 7/ cos a } 
we see that when = 0, a = B 7 , and when = 90°, a = B. Thus as the 
azimuth increases from one of these values to the other, the curve for a must 



1907.] On Light Elliptically Polarised by Reflexion. 



501 



pass from that representing B', which rises steadily as the angle of incidence 
increases, to that representing B, which is equal to B' at direct and grazing 
incidence, but falls to a minimum in the interval at the Principal Incidence. 

All these statements as to the variation of a are illustrated in the table and 
figures that follow. The table gives the values of a for various incidences 
and azimuths, the first entry referring to realgar, and the second to diamond. 
Figs. 14 and 15 below represent the same results graphically. 



e. 



£ = 10°. 



30°. 



60°. 



65°. 



67° 30'. 



67° 40'. 68°. 



70 c 



75°. 



80°. 



{ 
{ 
{ 
{ 

{ 
{ 



10 
30 

45 
60 



85 
89 



-4201 
-4218 

-4196 
-4194 

-4152 
-4169 

-4124 
-4141 

-4104 
-4128 

0-4098 
-4125 

-4098 
-4113 



-4613 
-4641 

-4421 
-4541 

-4171 
-4196 

-3912 
-3937 

-3712 
0*3737 

-3645 
-3717 

*3637 
'3659 



'6256 
-6328 

0*5558 
'5562 

-4552 
-4620 

-3576 
'3394 

'2012 
*2084 

-1378 
-1422 

-1292 
-1319 















s 












6761 
6671 

5921 
5932 

4842 
4871 

3546 
3453 

1845 
1834 

0987 
0779 

0918 
0527 



-6950 
-6984 

-6122 
-6156 

-5005 
-5018 

'3538 
-3551 

0-1832 
-1841 

-0789 
-0619 

-0792 
'0179 



-7004 
-6834 

-6144 
-6149 

-5018 
-5028 

*3550 
-3547 

-1837 
0*1841 

0-0796 
-0619 

-0798 
-0172 



-7028 
-6982 

-6170 
-6186 

-5042 
-5054 

-3565 
-3569 

-1848 
0-1849 

-0817 
-0627 

-0818 
-0197 



-7224 
-7310 

-6359 
0*6362 

0-5205 
-5212 

-3688 
-3703 

-1979 
-1962 

-1069 
-0809 

-0995 
-0528 



-7756 
-7823 

-6898 
*6926 

-5712 
-5753 

-4221 

0*4274 

-2668 
-2722 

-1968 
-1976 

-1875 
0-1865 



-8435 
0*8453 

-7619 
-7628 

-6561 
-6580 

-5313 
-5330^ 

-418& 
-4186 

-3756 
-3739 

-3667 
-3681 




Fig. '14.-TrEealgar. 
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*p 10 



£0 



30 40 50 

Fig. 15. — Diamond. 



60° , 70° 



80^ 



90" 



We have next to consider the mode of variation of the other semi- 
axis b, with the incidence and azimuth. The formulae already obtained 

give 

b 2 = a 2 tan 2 7 = (26)" 1 E 2 sin A sin 20 tan y = *' sin 20 tan y, 

where k is a constant for a given incidence. This enables us to discuss the 
variation of b with the azimuth. We see that 6 vanishes with y when = 
and = 90°, so that we should expect it to pass through a maximum some- 
where in the interval. The above equation for b 2 gives 

cos 2d" 



2hdh = 2 cos 20 tan 7 .d* + S i^. dy = ^23^ 



fC 



cos 2 7 



cos 2 7 



cos 2#~f- 



cos 2y. 

Thus 5 is a maximum when cos 20 = — cos 2a/cos 27. Except for the 
change of sign, this is the same equation as that dealt with above on p. 499. 
It corresponds, in the notation of that page, with the equation 

4:SX ^ * 



(1— x)(c+x) _ _ J 1 
(l+aO(c-a) L 



(1 + a?) 2 . 

In the present case we require 1 — ^/c — x to be negative, and, as before, we 
have 



■x __ — cos A 
— -j 



•# 



Hence below the Principle Incidence we must take the upper sign, and the 
lower sign beyond it. We thus obtain 

, 2 n 1 + e cos A 
tan 2 = _-—- _, 

e(e+ cos A) 
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the upper or lower sign being taken according as the incidence is less or 

greater than the Principle Incidence. In both cases we get- a real value 

for 8 for any given incidence. The magnitude of the corresponding 

maximum is 

6 = Esin A(l±2ecos A + e 2 )-i 

A simple geometrical construction for tan 2 # and the maximum of b is 
easily obtained. In fig. 16, take AB = 1, AC = e, and the angle BAO equal 




Fig. 16. 

to 7T — A or A, according as the incidence is lessor greater than the Principal 
Incidence. Draw AD perpendicular to BC, and in it take a point P so that 
AP = E' x AD. 

Then AP represents the maximum value of 6, and BD/CD represents tan 2 0. 
To prove this, draw BE and CF perpendicular to CA and AB respectively. 
Then 

sin A (1 + 26 cos A + e 2 )-* = BE/CE = sin C and E = e . B' = E' . AC. 
Hence 

Z>= EsinA(l + 2ecosA + e 2 )-£= K'.ACsinC = B'.AD; 

1 + e cos A BF 



also 



tan 2 6 = 



e(e±cosA) AC . CE ' 

but BF = BF . AB = BC . BD, and AC . CE = CD . CB, so that 

tan 2 6 = BD/CD. 

The changes in the magnitude of 6 and b are readily discussed by means of 
this construction. As A increases from to Jtt, BC/CD increases from 
unity. It afterwards diminishes to unity, so that 6 increases from 45° to XI, 
and then diminishes to 45°. Also, as A increases from to |tt, AD increases 
from zero to e (1 + e 2 )-*, and E' also increases, so that b increases steadily in 
this range from to E (1 + e d 2 )~K After this AD diminishes to zero. In 
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the neighbourhood of the Principal Incidence we have E sin A = K ,* so that 
in this region b = E /BC. As we go from the Principal Incidence, BC 
increases, so that b diminishes. Thus the maximum value of b is found at the 
Principal Incidence, and its value is E (1 + e 2 )~i The following table 
gives the position and magnitude of the maximum value of b for various 
incidences:- — 



(f>. 


10°. 


30°. 


60°. 


65°. ' 


67° 30' 


6 (Realgar) 


45° 22' 
0-0011 

45° 22' 
. 0-0002 


48° 30' 
-0107 

48° 29' 
-0023 


66° 19' 
-0521 

65° 40' 
0-0112 


73° 46' 
-0714 

75° 
0-0143 


83° 17' 
0-0786 

85° 33' 
-0170 


b (Realgar) 





4>. 


67° 40'. 


68°. 


70°. 


75°. 


80°. 


9 (Realgar) 


82° 48' 
-0785 

88° 37' 
0-0172 


80° 51' 
-0783 

84° 14' 
0-0172 


73° 40' 
-0769 

75° 34' 
-0161 


64° 40' 
-0625 

64° 11' 
-0155 


56° 47' 
-0539 

56° 45' 
3 -0119 


b (Realgar) 





The ridge of maxima is represented in plan and elevation in figs. 17 
and 18 below. 




icr ^(^ 3(r wr so w 60 zo 8tr 9(r 

Fig. 17. — The continuous curve refers to realgar, the dotted curve to diamond. 

* See ''Boy. Soc. Proc., 5 A, vol. 76, 1905, p. 61. 
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w 



70 



80 1 - 



-? io u eo° 30 p 40° 50° 

Fig. 18. — The continuous curve refers to realgar, the dotted curve to diamond. 



90 



In considering the variation of h with the incidence for a given azimuth, we 
take h 2 = (2e) _1 K 2 sin A . sin 20 . tan 7 = |E' E sin A., sin 20 . tan 7. 

As the incidence increases from direct incidence to the Principal Incidence, 
EsinA increases throughout, except at the end where it is stationary.* 
Tan 7 and E' also increase, so that h increases from zero upwards. When the 
Principal Incidence is passed, 7 decreases again to zero, so that h does likewise. 
It is evident then that h passes through a maximum between the Principal 
Incidence and grazing incidence, and that, as the azimuth increases to 90°, 
the position of this maximum draws closer and closer to the Principal 
Incidence. The following table gives the values of b for realgar and diamond 
for various incidences and azimuths : — 



0. 



4>=io°. 



30°. 



60°. 



65°. 



67° 30'. 



67° 40'. 



68°. 



70°. 



75°. 



80°. 



10° 
30° 
45° 
60° 

75° 
85° 
89° 



{ 

r 
1 

{ 

{ 
{ 
{ 
{ 



-0004 
-0001 

-0010 
0/0001 

o-oon 

-0002 

-0010 
-0001 

-0006 
-0001 

-0004 

-oooi 






-0032 
-0007 

-0086 
-0019 

-0106 
-0023 

0*0098 
-0021 

-0059 
0*0013 

0-0021 
-0009 

-0004 
-0001 



-0102 
0-0023 

-0299 
-0056 

-0417 
-0084 

-0515 
-0105 

-0475 
-0104 

0-0241 
-0053 

-0052 



-0134 
-0027 

-0378 
-0076 

-0517 
-0108 

-0632 
-0131 

-0665 
-0143 

-0462 
0-0117 

-0100 



0-0012 0-0035 



























0138 
0030 

0396 
0085 

0560 
0118 

0686 
0149 

0765 
0169 

0617 
0170 

0123 
0118 



0-0139 
-0030 

-0397 
0-0085 

-0562 
-0119 

-0690 
-0150 

0-0770 
-0168 

-0618 



-0139 
0-0031 

0-0402 
0-0086 

-0569 
-0120 

-0698 
-0150 

-0778 
-0168 

-0611 



0-0170 0-0170 



-0124 
-0124 



-0123 
0-0112 



*0143 
-0032 

-0411 
-0087 

-0585 
-0124 

-0737 
-0150 

-0766 
-0159 

-0495 
-0136 

'0107 
0-0042 



-0123 
-0034 

-0378 
-0095 

0-0527 
-0132 

-0617 
-0153 

0-0564 
0-0130 

-0266 
-0066 

-0056 
-0014 



-0082 
0-0026 

-0250 
0.-0082 

-0430 
0-0111 

-0480 
-0119 

-0334 

0-0087 

-0152 
-0034 

-0039 
0-0006 



* See 'Koy. Soc. Proc./ A, vol. 76, 1905, p. 61. 
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These results are represented in figs. 19 and 20 below. 

008 




+ 50 



60" 
Fig. 19. — Eealgar. 



70 v 



80 




Fig. 20. — Diamond. 



The area of the elliptic orbit is 

7T ab = (2e)~ l ttE 2 sin A sin 20 = JttB'R sin A . sin 20. 
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0*89 K 



Fia. 21. — Kealgar. 



Sb 
n 







k 



&*89 C 



6=9Cr 



Fig. 22.— Diamond. 
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For a given incidence this is proportional to sin 20, it therefore grows 
from zero to a maximum when 6 = 45°, and diminishes again to zero as 
6 approaches 90°. 

For a given azimuth the area is proportional to E'B sin A. As <p increases 
from zero to the Principal Incidence, E' increases, and E sin A increases, 
except near the Principal Incidence, where it is stationary. Hence the area 
of the ellipse constantly increases within this range. In the neighbourhood 
of the Principal Incidence E sin A is constant, so that the area is proportional 
to E', and continues to increase throughout a short range beyond the 
Principal Incidence. When this region is passed, sin A diminishes rapidly to 
zero, so that the area diminishes to zero, passing through a maximum between 
the Principal Incidence and grazing incidence. 

Figs. 21 and 22 represent half of the elliptic orbit, drawn to scale for 
various azimuths at the Principal Incidence, 



